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Abstract: We analyze the Euclidean version of supersymmetric quantum mechanics 
on the lattice by means of a numerical path integral. We consider two different lattice 
derivatives and improve the actions containing them with respect to supersymmetry by 
systematically adding interaction terms with non-zero extent. To quantize this improve- 
ment, we measure boson and fermion masses and Ward identities for the naive as well as 
the improved models. The masses are degenerate in all models, but the magnitude of the 
Ward identities decreases significantly for both derivative operators using the improved 
actions. This is a clear sign that the breaking of supersymmetry due to lattice artifacts is 
reduced. 
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1 Introduction 

Supersymmetry (SUSY) is assumed to play an important role in the search for a uni- 
fied theory which incorporates all known physical interactions. As a variety of interesting 
aspects of supersymmetric theories like super Yang-Mills are not accessible via perturba- 
tion theory, other approaches are required. Promising candidates are numerical methods 
like the computation of path integrals on space time lattices which have been applied to 
quantum chromodynamics with great success. There is, however, a fundamental problem 
that supersymmetric theories on discrete space-time suffer from, which is the failure of 
the Leibniz rule on the lattice [1, 2]. The naive ansatz for a lattice action usually is not 
invariant under infinitesimal translations, which it would have to be in order be fully su- 
persymmetric. Nevertheless, it is possible to formulate lattice theories which are not fully 
supersymmetric but describe the correct theory in the continuum limit. One way to ensure 
this is the orbifolding procedure [3] that allows to keep an exact subgroup of the original 
symmetry on the lattice. Reviews on the progress made with this method are given in 
Refs. [4, 5]. 
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In this work, we consider supersymmetric quantum mechanics (SUSYQM) on the lat- 
tice, which is a very simple model and thus a good playground to test new ideas. For 
SUSYQM, degenerate boson and fermion masses are obtained at non-zero lattice spacing 
and with interaction by including the Wilson mass also in the bosonic derivative operator 
[6]. The quadratic terms in the action are then exactly supersymmetric and the only ex- 
plicit breaking comes from the interaction term which connects bosons and fermions. This 
action also yields a continuum limit for the masses which is compatible with the exact 
value. Further improved actions for SUSYQM are discussed in Ref. [7], where a lattice 
action preserving one exact supersymmetry, introduced in Ref. [6], is numerically studied. 
Recently, SUSYQM on the lattice has been treated by means of a fermion loop formulation 
[8]. It is also possible to construct fully supersymmetric lattice theories, but this comes at 
the cost of a non-local action [9], see also [10, 11]. Besides, it is not clear how to generalize 
this construction to lattice theories which contain non-polynomial functions of fundamental 
fields like gauge links. Another approach is to keep only a remnant of the continuum SUSY 
on the lattice which is motivated by a blocking ansatz [12], but no lattice action has been 
found for this case so far. 

We will follow none of these approaches but yet another route, which is focused on 
finding a compromise between locality and SUSY. To this end, we allow the interaction 
terms to combine fields over a finite extent (still being ultra-local) and optimize the mag- 
nitude of these combinations to minimize the breaking of translational invariance. Due 
to the connection of SUSY and Poincare symmetry, we assume that this also leads to an 
improvement of supersymmetry, which we will confirm numerically. Contrary to most of 
the other approaches, we want to improve SUSY at finite lattice spacing rather than fo- 
cusing on the continuum limit. However, a combination of these improvements could be a 
promising ansatz for future lattice calculations. 

We organize our work as follows. In Sec. 2, we introduce the continuum and lattice 
versions of supersymmetric quantum mechanics. In Sec. 3, we discuss our method to 
construct lattice models of SUSYQM with finite extent interaction terms. Concrete lattice 
models are introduced in Sec. 4, which are classified by the two different lattice derivatives 
we consider. We present our numerical results in Sec. 5, which include boson and fermion 
masses that are degenerate for all the models. We furthermore show that the improved 
models have smaller Ward identities, which indicates a smaller breaking of SUSY. We 
summarize our findings and conclude in Sec. 6. 

2 Supersymmetric quantum mechanics 
2.1 Continuum 

This section contains a very brief introduction to SUSYQM in the continuum. The eu- 
clidean action of the theory reads 
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with bosonic real fields x an d F, and fermionic fields tp and tp which depend only on 
the real time variable t. dt denotes the derivative with respect to t and W{x) is a yet 
unspecified function called superpotential. The action S is invariant under supersymmetry 
transformations generated by 



M 



( 1 \ 
-d t 
-d t -10 
\ / 



and M 



/o -1 o\ 
-dt 


\d t -1 0/ 



(2.2) 



which have to be understood as linear operators acting on the field multiplet (x, F, ip, ip) T ■ 
The anti-commutator of these generators is proportional to the derivative dt, i.e. 



{M, M} = 2d t 



(2.3) 



As the auxiliary field F occurs at most quadratically in the action, it is readily integrated 
out, resulting in the substitution F — > —W(x), and consequently the on-shell action 



dt 



1 



dW 



1 



8x * + |W W 



(2.4) 



In the on-shell formulation, the variation of the fields under SUSY transformations gener- 
ated by M and M is given by 



M : # = 0, 

M : 5ip = edtx + eW . 



Sip 



-edtx + eW , 



0. 



5x = etp, 
$X = -tip , 



(2.5) 
(2.6) 



where e and e are infinitesimal Grassmanian parameters. To obtain a specific model of 
SUSYQM, we choose the superpotential 



W(x) =mx + gx 3 , 
with the mass m and the coupling constant g and obtain 



5° 



dt 



■^xdtX + ^dtip + mipip + Zgipx 2 ^ + ~ ( m X + #X 3 ) 2 



(2.7) 



(2i 



The theory resulting from this action has a supersymmetry which is not spontaneously 
broken, i.e. the first bosonic and fermionic excitations have the same energies. 



2.2 Lattice 

A lattice version of SUSYQM is obtained by discretizing the time variable t. We denote 
the number of lattice sites by n and the lattice spacing by a. All the fields acquire an index 
that denotes the lattice site and doubly occurring indices are implicitly summed from to 
n — 1. We consider the lattice version of the continuum action given in Eq. (2.1), 



1 



Xi(^ 2 )ijXj + i>Nij'ipj 



1 - dWi 

-FiFi + ^i '-tPj-FiWi 

2 Xj 



(2.9) 
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with the superpotential 

Wi = rhijXj + gT ijki XjXkXl ■ (2-10) 

The tensor T and the coupling constant g define the interaction term, V is a lattice deriva- 
tive operator and m is a general mass term. Each of the three will be specified later as we 
construct concrete lattice models. After integrating out the auxiliary field F, we obtain 
the lattice on-shell action 



5° 



\xi{^) 2 ijXj + i>iViji>j + ipimjipj + 3gT ijM ifjiXjXk^i + \ WiWi 



(2.11) 



The discretized SUSY generators M and M are obtained from their continuum counterparts 
defined in Eq. (2.2) via replacing the derivative dt by the lattice derivative V. In the on- 
shell formulation, they yield the variation of the fields 

M: Sipi = 0, 6%bi = -eVi jX j + eWi, 5 X i = &pi , (2.12) 

M: 6$i = eVijXj + eWi, 6^ = 0, S X i = ~&l>i ■ (2-13) 

The variation of the action (2.11) with respect to these generators is 

M : eSS = egipi(3T ijk iXjXk^lmXm -^ijTjkimXkXlXm) , (2.14) 
M : e5S = egipi(3Ti jk iXjXk^imXm -^ijTjkimXkXiXm) ■ (2.15) 

Contrary to the continuum theory, these variations do not vanish, as V does not obey a 
Leibniz rule in general. 



3 Method to construct a SUSY-improved lattice action 

Due to the connection of supersymmetry and Poincare symmetry, a fully supersymmetric 
action has to be invariant under arbitrary shifts of the fields generated by the lattice 
derivative V. We denote this property of a lattice action by continuous translational 
invariance (TI). This is an additional requirement which is not equivalent to discrete TI, 
i.e. the invariance of the action under shifts of one lattice spacing. The latter is a common 
property of generic lattice actions and is fulfilled by all actions considered in this work. As 
a consequence of discrete TI, quadratic terms in the action also fulfill continuous TI if the 
lattice derivative V is anti-symmetric, which is usually the case. However, for higher order 
interactions as characterized by the tensor T in the superpotential (2.10), continuous TI 
is one of the key problems of lattice SUSY. For the momentum space representations of 
the lattice derivative V and interaction tensor T, as defined in App. A, this invariance is 
equivalent to 

[V(p) + V(q) + V(r) + V(-p -q-r)] T(p, q,r)=0 (3.1) 

for all possible values of the momenta p, q and r. This equation cannot be fulfilled by a 
local interaction term [9]. Note that formally, Eq. (3.1) is solved by T(p,q,r) = 5\V(p) + 
V(<7) + V(r) + V(— p — q — r)] (here, 5 denotes the Kronecker symbol), as has been analyzed 
for V(p) ~ sin(ap/2) [1, 13]. However, on a finite lattice which is required for numerical 
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Table 1. Most local contributions t" bc to the interaction tensor T lmp and their locality L, Eq. (3.5), 
up to L = 4. 

simulations, such a choice of T usually does not yield a sensible interaction term, because 
the argument of the Kronecker symbol only becomes zero for very few combinations of the 
momenta p, q and r [1]. 

Our method will be to minimize the left hand side of Eq. (3.1) with a tensor T that 
is allowed to spread over a few lattice sites. This can be seen as a compromise between 
locality and supersymmetry. The most local interaction is the point-like one called 'naive', 

Tmaive sr /n ri\ 

ijkl = d ijkl , (3.2J 

with 5ijki = 5ij5ikSu. 

General contributions t to the interaction tensor T can be characterized by three lattice 
shifts a, b and c (due to discrete TI), and read 

(t abc ) ijkl = SYMMT {5 

i,j+a,k+b,l+c } , (3-3) 

where SYMMT{} means a symmetrization in the indices i, j, k and I, and under time 
reversal 1 . Due to these symmetrizations, the conditions 

0<a<5<c,and a<c-b, (3.4) 

can be imposed on the shifts without loss of generality. As a quantity to classify such 
contributions, we use their degree of locality L, which we define by the quadratic deviation 
from their center of mass X, 

L(t abc ) = (0-X) 2 + (a-X) 2 + (b-X) 2 + (c-X) 2 , X = (0 + a + b + c)/4. (3.5) 

In Tab. 1 we list the most local contributions up to L = 4, this list can easily be extended. 
The ordinal number of each tensor is denoted by an upper index a in t a . These seven 
tensors are exactly those, where fields are multiplied from lattice sites that are at most 
two sites apart, i.e. c < 2. We restrict ourselves to these tensors because we assume that 
they suffice to check our method to minimize the breaking of SUSY. However, one can in 

1 Time reversal means a substitution of {i,j, k, I) — > (— i, — j, —k, —I), where negative indices are projected 
back to the range (0, n — 1) via —i — > —i mod n. 
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principle allow for arbitrary contributions to the interaction tensor. In the simulations we 
will use an interaction tensor T imp that is the sum over the contributions t a , 

T™* = Y,c a t a , (3.6) 

a 

with (real) weights c a subject to the constraint 

J> = 1. (3.7) 



a 



The latter ensures that the superposition does not change the overall amplitude of the 
interaction tensor, which is given by the coupling constant g. 

As a quantity that measures the SUSY breaking of the interaction T (in combination 
with the lattice derivative V), we choose 

t/o 2 
B(T,V)= Yl \[V(p)+V(q) + V(r)+V(-p-q-r)]T(p,q,r) (3.8) 

p,q,r=— Tr/a 

i.e. the quadratic sum over all momenta on the left hand side of Eq. (3.1). We shall 
minimize this quantity with respect to the weights c a . Using a Lagrange multiplier A to 
ensure the constraint (3.7), this can actually be done analytically, because the resulting 
function 

B(T, V) + c a -l)=J2 Y aP c a cp + A c « " A ( 3 - 9 ) 

a a, (3 a 

is linear and bilinear in (c Q ,A) and thus the conditions for extrema are a set of linear 
equations for them. In each of the cases we consider, these have a unique solution for the 
c a . The involved coefficients 

t/o 2 
Y aP = Yl \[V(p)+V(q)+V(r) + V(-p-q-r)]\t a (p,q,r)tP(p,q,r) (3.10) 

p,q,r=— Tr/a 

are real (because the tensors t a are real in momentum space), and fulfill Y a p = Y@ a . They 
depend weakly on the number of lattice sites n. In the following, we therefore evaluate 
them in the thermodynamic n — > oo limit, where the sums over momenta are appropriately 
replaced by integrals. Taking the continuum limit a — > actually makes no difference 
because the lattice spacing can be factorized out of B for the lattice derivatives we consider, 
so it does not play a role in the minimization. 

Because of the symmetrization, the functions t a (p,q,r) consist of combinations of 
cosines, e.g., 

t 2 (p, q, r) oc cos(ap) + cos(aq) + cos(ar) + cos(ap + aq + ar) . (3-11) 

The antisymmetric lattice derivatives are odd functions in the momenta. For the symmetric 
derivative to be used in Sec. 4.2, V(p) is a sine function and the integrals occurring in Y a p 
can be calculated analytically. 
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As a rough estimation for the achieved improvement of the interaction tensor, we 
introduce the (non-physical) quantity 

Q^) = B ^^y ( 3 - 12 ) 

where T imp is the interaction tensor that results from a minimization of B. Values of Q for 
the improved models are provided in the respective sections. 

4 Lattice discretizations 

We consider four different discretizations based on the action given in Eq. (2.11). These 
contain different kinds of lattice derivatives and interaction terms which are described in 
detail below. 

4.1 Naive and improved actions with the SLAC-derivative 

The derivative used in the first models is the SLAC derivative, which is defined in momen- 
tum space by 

V SL Ac(p) = ip. (4.1) 

This choice makes a Wilson mass obsolete as the SLAC operator has no doublers, which 
however comes at the cost of being non-local. Therefore, the mass matrix becomes trivial, 

rhij = mbij , (4.2) 

with the bare mass m. 

At first, we introduce the 'naive SLAC model', whose interaction tensor is given by 
T naive of Eq. (3.2). Secondly we construct an improved model with an interaction tensor 
that contains contributions up to L = 4, or equivalently, with products of fields which are 
up to two lattice sites apart, 

7 

T im P = £ Cq t a ( 43 ) 
a=l 

The coefficients c a were determined as described in Sec. 3, resulting in 

ci = 0.0754; c 2 = 0.3389; c 3 = 0.2057; c 4 = 0.1687; 

c 5 = 0.1597; c 6 = 0.0421; c 7 = 0.0095 . (4.4) 

We refer to this choice of the interaction tensor as the 'improved SLAC model'. This model 
yields an improvement of 

Q(V SLAC ) « 2.3 • 10- 6 (4.5) 

as defined in Eq. (3.12). 
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4.2 Naive and improved actions with a Wilson mass 

In these models, we insert a symmetric difference operator into the action, which is defined 
by 

l 

(V symm )ij = — (<S,-,i+i - <5j,i-i) • (4.6) 
A Wilson mass has to be included to remove the doublers, so the mass matrix is 

TCL 

rhij = m5ij - — - 2<% + , (4.7) 

where we chose r = 1 for the Wilson mass parameter. For the 'naive Wilson model', 
we insert the naive interaction term as defined in Eq. (3.2). The contributions to the 
improved interaction term are again the ones from Eq. (4.3). This time, the constants that 
are obtained by a minimization of B(T, V symm ) are 

d = 79/1907,; c 2 = 448/1907; c 3 = 318/1907; c 4 = 432/1907; 

c 5 = 384/1907 ; c 6 = 216/1907 ; c 7 = 30/1907 ; (4.8) 

We call the model that results from this choice the 'improved Wilson model'. Here, we 
have obtained an improvement of 

Q(V syimn ) « 5.2 • 10" 4 . (4.9) 

5 Numerical results 

In this section, we give the numerical results for the masses and Ward identities we have 
measured with the actions defined in the previous section. For each model we have chosen a 
bare mass of m = 10 and a volume of na = 1, which ensures that the Compton wavelength 
of the lowest mode fits easily into the volume. All the field configurations were obtained 
by a hybrid Monte Carlo algorithm similar to the one used in Ref. [9]. 

5.1 Masses 

According to previous works, the naive SLAC model [14] as well as the naive Wilson model 
[6] defined in Sec. 4 should have degenerate boson and fermion masses. In order to check 
this and to see if it holds also for the improved models, we have calculated these masses 
for all the models we have defined. This has been done by linear fits to the respective 
propagators in a logarithmic representation as described in Ref. [9]. Unlike there, we have 
used a simple Gaussian filter 2 to smoothen the fermionic propagator for the SLAC models, 
because more complicated filters have given only negligibly different results. We have 
chosen a coupling constant of g = 100, which is in the regime of strong coupling, to be 
able to compare our results with the previous works mentioned above, which have used the 
same values of g and m. For each lattice size and model an ensemble of 10 6 configurations 
has been created, which is sufficient to obtain small errors for the propagators. As can be 

2 Defined by the replacement of the propagator (^oV'O — > c 12i e _( - fe_i ' ) ^ 2 (tpotpi) , with c = 1/X^fc e_fe ^ ■ 
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Figure 1. Masses for the SLAC models (left) and Wilson models (right) for m — 10 and g = 100. 
The linear fits were obtained by minimized square deviation. The horizontal dashed line is the 
exact continuum mass (= 16.865). 
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continuum 
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continuum 
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naive SLAC 


16.784 (94) 


16.870(5) 


improved SLAC 


16.813(50) 


16.793 (4) 


naive Wilson 


16.815(50) 


16.778(5) 


improved Wilson 


16.650 (43) 


16.645 (4) 


exact 


16.865 


16.865 



Table 2. Continuum mass extrapolations for the various actions. 



seen in Fig. 1, fermion and boson masses are fairly equal inside the error bars for all models 
and lattice spacings. For the naive SLAC model, the masses are compatible to the ones 
obtained in Ref. [14], where the same model has been considered. This holds also for the 
masses of the naive Wilson model, which have been calculated in Ref. [6]. The values of 
the respective continuum interpolations, which have been determined by a linear fit with 
minimized square deviation, are given in Tab. 2. 

The naive models are closer to the exact continuum mass 3 for all values of the lattice 
spacing, the reason for which is presumably the larger extension of the interaction term in 
the improved models, as the extension is zero in the continuum. This seems to be the price 
one has to pay for the improved supersymmetry. However, the masses of the improved 
models also have a clear tendency towards the correct continuum mass, even if the linear 
fits do not hit the exact value in all cases. We assume that combining our method with 
additional improvements of the action like those introduced in [7] would result in a much 
faster convergence of the masses to the continuum limit. 

3 The (degenerate) exact continuum mass can easily be obtained by a numerical treatment of the Hamil- 
tonian of SUSYQM [14]. 
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Figure 2. Functions I\^{k) denned in Eq. (5.1) and (5.2) associated with the Ward identities for 
both SLAC models (left) and the improved SLAC model alone (right) for n = 21, m = 10 and 
g = 800. Note the different axis scaling. 

5.2 Ward identities 

For all models, we have computed expectation values belonging to the two Ward identities 
(Wis) derived in App. B, 

WI 1: h(k) = (4^o + (-VoiXj + Wo) X k) = (SSiPaXk) , (5-1) 
WI 2: I 2 {k) = + (VojJC, + Wo) Xk) = (SS^oXk) , (5.2) 

where we have set I = in Eqs. (B.7) and (B.8) without loss of generality, because these 
expectation values are only functions of the difference I — k due to discrete TI. They vanish 
if the action is exactly supersymmetric, i.e. 5S = 5S = 0, so their magnitude is a measure 
for the breaking of SUSY. In Figs. 2 and 3 we only show the expectation values related 
to 5S (5S), i.e. the right hand sides of Eqs. (5.1) and (5.2), because these suffer from 
far smaller numerical errors than the other expectation values. We have checked for each 
Ward identity that both expectation values coincide within the error bars. Due to the 
connection between the Ward identities mentioned in App. B, these expectation values 
fulfill I\{k) = l2(—k). We have chosen a very strong coupling of g = 800 and a small 
lattice size of n = 21 for all the models to have a strongly broken SUSY in the naive 
models. 

As shown in Fig. 2, the scale of the Ward identities in the improved SLAC model 
is smaller by four orders of magnitude compared to the naive model. The number of 
configurations is 10 for the naive SLAC model and 10 8 for the improved SLAC model. In 
the latter case, better statistics have been required to get reliable results due to the small 
value of the Ward identities. 

The same Ward identities have been computed for the Wilson models. These are 
shown in Fig. 3, where again 10 7 field configurations have been evaluated for the naive 
model and 10 8 for the improved model. There clearly is an improvement of about one 
order of magnitude, which however is much less than in the SLAC model. One possible 
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Figure 3. Functions h,2{k) denned in Eq. (5.1) and (5.2) associated with the Ward identities for 
both Wilson models (left) and the improved Wilson model alone (right) for n — 21, m — 10 and 
g = 800. Note the different axis scaling. 

reason for this fact are the different values of Q that were achieved, because these differ by 
two orders of magnitude. Presumably, the SLAC derivative is also more suited to allow for 
a supersymmetric interaction term in the first place. This assumption is backed by the fact 
that the only possible construction of an exactly supersymmetric (albeit non-local) lattice 
action with interaction requires the use of the SLAC derivative [9]. 

6 Summary 

We have constructed lattice actions of supersymmetric quantum mechanics which are im- 
proved with respect to SUSY and have compared them to naive discretizations. One type 
of models contains the SLAC derivative in the kinetic terms of the action, while the other 
incorporates the naive symmetric difference operator and consequently a Wilson mass term. 
We have made a compromise between supersymmetry and locality by allowing the interac- 
tion term of the improved actions to connect fields within two lattice spacings or less. Via 
numerical lattice simulations we have determined boson and fermion masses which have 
turned out to be degenerate for each model. The Ward identities are much smaller for the 
improved models than for the naive ones. We therefore conclude that the improved actions 
indeed have better properties with respect to supersymmetry. 

The method we have developed to construct improved lattice actions can easily be 
generalized to include interaction terms that extend over more than two lattice sites or 
are polynomials of higher rank. Furthermore, our method could also be applied to higher 
dimensional Wess-Zumino type models. This would not render numerical simulations too 
expensive, because the interaction terms we have used are ultra-local. Therefore, our 
technique could be an ingredient for the construction of supersymmetric lattice actions 
of more realistic models than SUSYQM, possibly in combination with other improvement 
methods. 
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A Fourier transforms 

We define the Fourier transform of a general tensor G with a number of M indices by 

n-l 

G(pi, ..,p M ) = n-f Yl G h,~,iM e wE2ilftJs • (A.l) 
iivjM=o 
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In this work, we consider only tensors which fulfill discrete TI, i.e. they are invariant under 
an increase of all their indices by one, in which case we obtain 

/ M \ n-1 

G(pi, .., PM ) = S P I J> J n-f +1 Yl G h,-dM-uO e ia ^^ , ( A .2) 

\»=1 / jivjM-i=0 
where <5p is a periodic <5- function, defined by 

<5 P (g) = S (q mod ^ . (A.3) 

As the 5 p-function constrains the sum of the momenta, we can drop one of them to arrive 
at the final form of the tensor G in momentum space, 

n-1 

G( Pl , .., PM -i) = n-f + 1 J2 G h,-,3M-^ e-^ 1 ^ . (A.4) 

jl,--jM-l=0 

It is implicitly understood that the omitted momentum pm is the negative sum of all the 
other momenta (projected to the first Brillouin zone). 

B Ward identities 

To obtain Ward identities, we consider the generating functional of a general lattice theory, 
defined by 

Z[J] = -1 J Dcf) e S[4>]+Jacfa f ( B-1 j 

with the action S, the fields <p a and the currents J a . a is a super-index which labels both 
field species as well as lattice sites. Doubly occurring super-indices are summed over all 
possible values they can take. N is the partition function defined by 

M = j Dcp e" 5M . (B.2) 

An infinitesimal transformation of the fields 

<Pa^<f>a + (6<f>) a = 0a + G ab (f> b (B.3) 

with the generator G yields 

Z[J] _». _L J Dcf) det [1 + G] e -S[<P]-SS[<p]+J a (<t>a+G ab 4> b ) 

« Z[J\{\ +tr[G]) + ^ J Dcf>(J a G ab cf> b - SSffl) e~ s ^ +J ^ , (B.4) 

where SS is the variation of the action with respect to the transformation. This has to be 
equal to Z[J], because the transformation is equivalent to a mere change of variables in 
the path integral. For a traceless G, we can therefore conclude 

J D4> (SS[<i>] + J a G ab fo) e ~ s ^ +J ^ = . (B.5) 
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Ward identities are obtained by deriving this expression with respect to currents and setting 
these to zero, 



d d 
dJ c dJ d 



D<t> (-5S&] + J a G ab <p b ) e-m+Je^ 

J D(j)(- 5S[<p](j) c <p d + G ca (j)a<Pd + 4>cG da <Pa) 

D<P ( - 6S[<t>]<l> c <t> d + {5<j>) c <j> d + MW)d) e" 5M = . (B.6) 



J=0 



In the case of SUSYQM, the generator G can be either eM or eM, whose actions are 
defined in Eqs. (2.12) and (2.13). We specify (f> c = Xk and (fid = ipi or tpi, where I and k 
label lattice sites, to obtain the Ward identities 

M: (^i + (-VijX j + W l ) Xk ) = (SSij l Xk), (B.7) 
M: (Mk + {VijX3 + W l )xk) = (6Sj>iXk), (B.8) 

where 5S and 5S are given in Eq. (2.14) and Eq. (2.15), respectively, and the average is 
defined by the path integral, i.e. 

There exists a connection between the two Ward identities, which can be shown as follows. 
We define 

W« = {Mi + (-V ljXj + Wi) Xk), (B.10) 
Wki = + {VijXj + Wi) Xk) ■ (B.ll) 

Due to the invariance of the lattice action S under a time-reversal of the x-field, i.e. 
a replacement x% ~^ X—i an d lattice translational invariance, it is easily checked that 
W k i = Wi k . 
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